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(54) Title: DETERMINING A PARAMETER OF A COMPONENT IN A COMPOSITION 
(57) Abstract 



A method of determining a parameter selected from the electrical conductivity and the volume fraction of a component in a composition 
comprising a plurality of components is provided. The method comprises measuring the electrical conductivity of the composition, and 
selecting a relationship between the conductivity of the composition and a plurality of composition parameters including, for each component, 
physical parameters representing the conductivity and the volume fraction of the component, said relationship being such that the components 
are substantially equally represented in said relationship by means of said physical parameters. The selected parameter of the component 
in the composition is determined from said relationship and the measured conductivity of the composition. 
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DETERMINING A PARAMETER OF A 
COMPONENT IN A COMPOSITION 

The present invention relates to a method of 
determining a parameter selected from the electrical 
conductivity and the volume fraction of a component in a 
composition comprising a plurality of components. The 
invention is of particular interest for determining the 
volume fraction of a component of an earth formation, for 
example to determine the hydrocarbon-content of a 
hydrocarbon bearing earth formation. Various known 
methods of determining such content have been applied 
until. now, which known methods are generally based on 
empirical models. 

One such known method is described in "Electrical 
conductivities in oil-bearing shaly sands", Waxman M.H. 
and Smits L.J.M., SPE paper 1863-A presented at 42nd Ann. 
Fall Meeting, Houston, October 1-4, 1967. 

This publication discloses a method of determining a 
parameter selected from the electrical conductivity and 
the volume fraction of a component in a composition 
comprising a plurality of components, wherein the 
electrical conductivity of the composition is measured, 
and a relationship between the conductivity of the 
composition and the conductivity of a component is 
selected. 

This known method uses the following relationship 
which is generally referred to as the Waxman-Smits model: 

C Q = C w / F* + BQ V / F* 

where C Q = conductivity of fully brine saturated rock 

C w = conductivity of brine present in the formation 
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F* = a formation factor 

B = equivalent conductance of sodium clay-exchange 

cations as a function of C w 
Q v = cation exchange capacity per unit pore volume . 
5 The results achieved with this known method are not 

always sufficiently accurate, probably because of the 
empirical nature of the Waxman-Smits model which provides 
a relationship between the earth conductivity and the 
various other parameters. 
10 It is, an object of the invention to provide a more 

accurate method of determining a parameter selected from 
the electrical conductivity and the volume fraction of a 
component in a composition comprising a plurality of 
components . 

15 The method according to the invention thereto 

comprises : 

measuring the electrical conductivity of the 
composition; 

selecting a relationship between the electrical 

20 conductivity of the composition and a plurality of 

composition parameters including, for each component, 
physical parameters representing the electrical 
conductivity and the volume fraction of the component, 
the components being substantially equally represented in 

25 said relationship by means of said physical parameters; 

and 

determining said selected parameter from said 
relationship and the measured conductivity of the 
composition . 

30 It is to be understood that by the electrical 

conductivity is meant the electrical conductivity itself 
or any quantity derived therefrom, such as the electrical 
resistivity. Furthermore, the feature that the components 
are substantially equally represented in the relationship 

35 implies that each component is represented in the 
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relationship in substantially the same way as any other 
of the components. 

With the method according to the invention it is 
achieved that results of increased accuracy are provided. 
5 The selected relationship takes accurately account of the 

individual contributions of the components to the 
conductivity of the composition. The relationship applied 
in the method according to the invention is symmetrical 
in the components, i.e. no component is favoured over 

10 another component. Moreover, it was found that the method 

according to the invention provides the desired accuracy 
for any percolation threshold of the components. In this 
respect it is to be understood that the amount of 
percolation of a component refers to the degree of 

15 continuity of the component in the composition. For 

example, vanishing percolation of a component implies 
that the component is fully dispersed in the composition, 
and full percolation of a component implies that the 
component is continuous throughout the composition. 

20 Advantageously said plurality of composition 

parameters includes at least one fitting parameter, and 
wherein each fitting parameter is determined by applying 
said relationship to a data set obtained by measuring the 
electrical conductivity of at least one sample 

25 representative for said composition for various 

magnitudes of at least one of said parameters. 

Preferably said plurality of parameters includes an 
auxiliary parameter depending on the geometrical 
configuration of the components in the composition. 

30 Accurate geometrical representation by the auxiliary 

component is achieved if said auxiliary parameter is 
selected so as to be a function of a plurality of 
variables, each variable depending on the conductivity of 
one of said components and a mixing coefficient, whereby 



BNSDOC3D: <WO 98296 16A1J_> 



WO 96/29616 PCIYEP96/01234 

ft 

- 4 - 

the mixing coefficients depend on the geometrical 
configuration of the components in the composition. 

Advantageously the step of determining each fitting 
parameter by applying said relationship to the data set 
5 of the component is carried out through an iterative 

process. Suitably the iterative process includes 
repeatedly applying said relationship in a minimisation 
scheme. The minimisation scheme is preferably applied to 
an incoherence between the measured electrical 

10 conductivities of said components and the electrical 

conductivities of the components as determined through 
said relationship . 

The invention will be described hereinafter in more 
detail and by way of the following example and 

15 comparative example. 

Example 

Consider an isotropic system with essential spherical 
inclusions in the form of an earth formation which 
essentially consists of four components: non-conducting 

20 porous rock matrix, non-conducting hydrocarbon fluid, 

conducting clay, and conducting brine. The conductivity 
of the formation depends on the fractional brine 
saturation of the pore space, and the hydrocarbon fluid 
component is grouped with the rock matrix, both being 

25 non-conducting. Thus, the hydrocarbon component and the 

rock matrix component only enter the equations with the 
sum of their volume fractions. The effective conductivity 
°eff of this earth formation is evaluated through the 
expression 

< G eff * °0> • (Lo e f f + (l-L)ao)" 1 = I ^ <°k " oq) . 
(La k + (I-Doq)" 1 



30 wherein oq represents the auxiliary parameter in the 

form of a conductivity tensor 
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k = 1 . . . N, N being the number of components 

a eff represents the conductivity tensor of the sample 

represents the conductivity tensor of component k 
<J>k represents the volume fraction of component k 
5 L represents the depolarisation tensor (shape 

tensor) 

Preferably the depolarisation tensor is positive and 
has unit trace* In an attractive embodiment the 
depolarisation tensor equals 1/3 times the unit tensor, 

10 The term oq denotes an auxiliary parameter which can 

be thought of as being an additional host medium into 
which components are added until the host medium has been 
completely replaced by the components so that no volume 
fraction is associated with the host medium. The 

15 existence of the host medium enables the model to be 

symmetrical in all its constituents: none of the 
components rock, clay or brine in the model is favoured 
over any of the other components. The dependence of oq on 
various parameters, yet to be determined, governs the 

20 percolation behaviour of the model. Setting oq = <*brine 

leads to the known Average T-matrix Approximation, also 
referred to as the generalised Clausius-Mossotti 
equation. This model has a clear asymmetry between the 
brine component and the other components since only the 

25 brine component will percolate, irrespective of its 

volume fraction. Selecting a self-consistent host medium 
conductivity, oq = o e ff, leads to the known Coherent 
Potential Approximation, also referred to as the 
generalised Bruggeman equation. This model is symmetrical 

30 in all components but has the drawback of requiring 

unrealistically high percolation thresholds for each 
component . 

In a suitable embodiment, the auxiliary parameter oq 
is selected as follows: 
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oq = Z ; for k = 1, 2, 3 

wherein hj,. represents the mixing coefficient tensor 
pertaining to component k, which tensor contains mixing 

.5 coefficients representing geometrical information on the 

spatial distribution of the components in the formation. 
These coefficients determine the connectivity, i.e. the 
amount of percolation of the individual components. The 
coefficients are non-negative and fulfil the 

10 normalisation condition: 

Z h k = 1; for k = 1, 2, 3 

The normalisation relation ensures that the resulting 
15 effective conductivity a e ff satisfies the Hashin - 

Shtrikman bounds, which are well known to those skilled 
in the art. 

Furthermore, a component with a vanishingly small volume 
fraction can not percolate, hence the corresponding 
20 connectivity parameter should vanish: 

lim hfc = 0; for 4)^ 0 



Suitably the mixing coefficient tensor is selected to 

be 

hk - *k4>k v k ( ZMfcV" 1 

wherein k,n = 1 ... N, N being the number of 
components in said plurality of components 
25 Xfc represents the percolation rate tensor pertaining 

to component k 

(J)^ represents the volume fraction of component k 
v represents the percolation exponent pertaining to 
component k 
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Suitably at least one of h^, and v forms a fitting 
parameter. 

A data set on 27 shaly-sand core samples has been 
used to test the invention, which data set is described 
5 in the above indicated SPE paper. This publication 

provides C G - C w curves on the core samples ranging from 
almost clean sand (Q v = 0.017 eq/1) to extremely shaly 
sand (Q v =1.47 eq/1). The samples contained Kaolinite, 
Montmorillonite and Illite, either in combination or 

10 separately in each sample. The characteristic 

petrophysical data of each sample are listed in the 
appended Table, in which 4> denotes the porosity of the 
sample, . k denotes the permeability of the sample, and Q v 
denotes the cation exchange capacity per unit pore volume 

15 of the sample. The conductivity of each sample in fully 

brine saturated condition was measured for eight to ten 
salinities of the brine. Furthermore, concentration 
membrane potential measurements were made of the samples. 
The parameters in this model were selected as 

20 follows: 

1) Brine; 

The volume fraction of brine, <J>t>, is determined by the 
porosity, the amount of clay-bound water, and the water 
saturation S w . The brine conductivity Ofc> (=C W ) is 
25 determined by the brine salinity and the brine 

temperature. The two percolation parameters, 7^ and v, 
are free parameters. 

2) Rock / Hydrocarbon; 

The volume of hydrocarbons, <|>hc' is determined by the 
30 total porosity, the amount of clay-bound water, and the 

hydrocarbon saturation 1-S W , while the volume of the rock 
matrix, 4> r , is calculated using the sum rule and the 
volume fractions. Both the rock and the hydrocarbon have 
vanishing conductivity. The percolation parameters X r and 
35 Xfr c of both components was set at value 1. The mixing 
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coefficient pertaining to rock/hydrocarbon h r / hc follows 
from the condition Z h^ = 1 . 
3) Clay; 

The volume of clay 4> c and the clay conductivity o c are 
5 free fitting parameters. The percolation rate X c was set 

at a value 0, which is a suitable choice for non- 
laminated clays- It furthermore appeared that an 
additional free parameter did not give a significant 
improvement of the model fit to the data set. 

10 The C Q - C w measurements were made for an extreme 

salinity range, namely a brine salinity between 1 - 300 
g/1. For a given sample the brine volume fraction varied 
only slightly over the whole salinity range. In view 
thereof the percolation parameter v was set equal to 

15 . unity in the test, thereby reducing the percolation 

parameter h^ to a constant, and reducing the number of 
free parameters' to three. 

For each sample, a fit to the C G - C w curve was made 
by minimising the relative incoherence defined as: 

20 



25 



I A 2 C Q - C w = I 

salinities \ C 



Co* calc Co, meas 



o, meas 



wherein 

c o,calc = the calculated conductivity of the fully brine 
saturated rock samples; 

c o,meas = the measured conductivity of the fully brine 
30 saturated rock samples; 

X = summation over the salinities. 

The results for the three fitting parameters 4> c , a c 
and hfc, and the relative incoherence are given in the 
appended Table. 
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10 



15 



20 



25 



Furthermore, the Table gives the incoherence between 
the membrane potential (vj/ ca i C ) determined by the method 
of the invention and the measured membrane potential (y 

meas) : 



The membrane potential is a particularly interesting 
quantity for being a direct, non-destructive, measure of 
the clay contribution to the overall conductivity, which 
has not been used to determine the fitting parameters. 

To illustrate the invention more specifically, 
reference is made to the following comparative example. 
Comparative example 

As stated above Ref . 1 discloses, apart from the data set 
on the 27 core samples, furthermore an empirical model 
which is generally referred to as the Waxman-Smits model. 
To compare the method according to the invention with the 
Waxman-Smits model, the relative incoherence between the 
measured conductivities and the conductivities found from 
the Waxman-Smits model, and the relative incoherence 
between the measured concentration membrane potentials 
and the concentration membrane potentials found from the 
Waxman-Smits model, were determined. These relative 
incoherences for all 27 samples are listed in the Table. 
In applying the Waxman-Smits model, use has been made of 
the well known expression: 

C Q = C w / F* + BQ V / F* 
with F* = <\>~ m 



I A2 Mp = £ 




Vcalc Vmeas 



Vmeas 
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where m is a free parameter (also referred to as the 
cementation exponent) , Q v is determined from sample 
measurements, just as the porosity <t>, and the standard B- 
chart has been used to calculate the salinity and 
5 temperature effects on the conductivity measurements. 

From a comparison between the incoherence values 
found by using the method according to the invention, and 
the incoherence values found by using the Waxman-Smits 
model, it is clear that the method according to the 

10 invention provides improved results . Especially the 

extremely low incoherence values for the concentration 
membrane potential, which values are moreover fairly 
constant over the entire Q v range, indicates that the 
method according to the invention provides results of 

15 increased accuracy. 

The method according to the invention can suitably be 
applied to determine the volume fraction of brine or 
hydrocarbon in an earth formation, whereby a well-log 
representing the electrical conductivity of the formation 

20 is provided. Such application can, for example, be 

carried out in the following manner. The well-log of the 
electrical conductivity of the earth formation is made 
using a logging tool lowered in a borehole formed in the 
earth formation. For an isotropic formation with 

25 components brine (subscript B) ,' clay (subscript C) , and 

non-conducting rock + hydrocarbon (subscript R/HC) the 
rock and the hydrocarbon are grouped together because of 
their vanishing conductivities. The selected relationship 
then is: 

a eff~ a o 3 o k - o Q 

— = X <j>k- 

a eff + 2 °o k=l <*k + 2a o 

30 wherein 
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oo = H h k o k 

h k - M>k v k ( SX n (J) n v n )-l 

in which 

gq represents the auxiliary parameter 
k, n = 1 . . . N, N being the number components 
G e ff represents the conductivity of the earth formation 
5 a k represents the conductivity of component k 

4 k represents the volume fraction of component k 
h k represents the mixing coefficient pertaining to 
component k; 

X k represents the percolation coefficient pertaining to 

10 component k 

v k,n represents the percolation exponent pertaining to 
component k, n 
Each component k has four parameters: 4)^, o k , X k , and 
v k , of which <t> B , a B and o r /hc are directly measured- 

15 Furthermore, X c = 0 for dispersed clay. From the sum 

rules h^/HC and 4>R/HC follow. Parameters which are yet to 
be determined are o c , X B , v B and 4>c- These parameters are 
determined through forward modelling on experimental 
data. Gq, X b anc j v B are invariable over the geological 

20 formation , while <J> B will be depth dependent. The 

experimental data for the parameter determination consist 
of well-log measurements from a brine containing zone, 
laboratory Formation Resistivity Factor (FRF) 
measurements and brine saturation experiments. The log 

25 information from the brine containing zone is used to 

correlate the local parameter 4>c to suitable logs/log 
combinations, as is known to those skilled in the art of 
well logging. oq, X b and <j> B and the correlation of q>c to 
suitable logs/log combinations can be used in hydrocarbon 

30 bearing formations. From the well-log, the above 

relationship and the indicated parameters, the volume 
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faction of brine and thus also the volume fraction of 
hydrocarbon is determined as a function of depth* 
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CLAIMS 

1. A method of determining a parameter selected from the 
electrical conductivity and the volume fraction of a 
component in a composition comprising a plurality of 
components, the method comprising: 

5 - measuring the electrical conductivity of the 

composition; 

selecting a relationship between the electrical 
conductivity of the composition and a plurality of 
composition parameters including, for each component, 
10 physical parameters representing the electrical 

conductivity and the volume fraction of the component, 
the components being substantially equally represented in 
said relationship by means of said physical parameters; 
and 

15 - determining said selected parameter from said 

relationship and the measured conductivity of the 
composition . 

2. The method of claim 1, wherein said plurality of 
composition parameters includes at least one fitting 

20 parameter, and wherein each fitting parameter is 

determined by applying said relationship to a data set 
obtained by measuring the electrical conductivity of at 
least one sample representative for said composition for 
various magnitudes of at least one of said parameters. 

25 3. The method of claim 1 or 2, wherein said plurality of 

parameters includes an auxiliary parameter depending on 
the geometrical configuration of the components in the 
composition. 

4. The method of claim 3, wherein said relationship is 
30 selected to be 
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< a eff " °0) ■ (La eff + (l-L)ao)" 1 = I (°k " a 0) • < La k + 
(I-Doq)" 1 

wherein oq represents the auxiliary parameter in the 
form of a conductivity tensor 

k = 1 . . . N, N being the number components 
°eff represents the conductivity tensor of the 
5 sample 

o k represents the conductivity tensor of component 
k 

<J> k represents the volume fraction of component k 
L represents a depolarisation tensor 
10 5. The method of claim 3 or 4, wherein said auxiliary 

parameter is selected so as to be a function of a 
plurality of variables, each variable depending on the 
conductivity of one of said components and a mixing 
coefficient, whereby the mixing coefficients depend on 
15 the geometrical configuration of the components in the 

composition . 

6. The method of any of claims 3-5, wherein said 
auxiliary parameter is selected to be 

o 0 = Z h k a k 

20 wherein oq represents the auxiliary parameter in the 

form of a conductivity tensor representative for the 
conductivity in the three principal directions; 

k = 1 ... N, N being the number of components; 
<*k represents the conductivity tensor of component 
25 k; 

h k represents the mixing coefficient tensor 
pertaining to component k; 

7. The method of claim 5 or 6, wherein said mixing 
coefficients are selected so that the sum of the mixing 

30 coefficients substantially equals unity. 

8. The method of any of claims 5-7, wherein said mixing 
coefficients are non-negative. 
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9. The method of any of claims 5-8, wherein each mixing 
coefficient is selected to be a function of at least the 
volume fraction of the component pertaining to said 
mixing coefficient. 
5 10. The method of claim 9, wherein said function is a 

monotonous increasing function in the volume fraction of 
the component pertaining to the mixing coefficient. 

11. The method of claim 9 or 10, wherein said function is 
selected so that the mixing coefficient vanishes for 

10 vanishing volume fraction of the component pertaining to 

. the mixing coefficient. 

12. The method of any of claims 5-11/ wherein each mixing 
coefficient is selected as 

h k = *-k<t>k v k < ^n^n)" 1 
15 wherein k, n = 1 . . . N, N being the number components 

in said plurality of components 

represents a percolation rate tensor pertaining to 

component k 

4*^ represents the volume fraction of component k 
2° v k,n represents a percolation exponent pertaining to 

component k, n 

13. The method of claim 12, wherein at least one of h^, X 
^ and v forms a fitting parameter. 

14. The method of any of claims 4-13, wherein the 
25 depolarisation tensor is positive. 

15. The method of any of claims 4-14, wherein the 
depolarisation tensor has unit trace. 

16. The method of any of claims 4-15, wherein the 
depolarisation tensor equals 1/3 times the unit tensor. 

30 17. The method of any of claims 1-16, wherein the step of 

determining each fitting parameter by applying said 
relationship to the data set is carried out through an 
iterative process . 
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18. The method of claim 17, wherein the iterative process 
includes repeatedly applying said relationship in a 
minimisation scheme. 

19. The method of claim 18, wherein the minimisation 
5 scheme is applied to a mismatch between the measured 

electrical conductivities of said components and the 
electrical conductivities of the components as determined 
through said relationship. 

20. The method of any of claims 1-19, wherein said 
10 composition includes an earth formation, 

21. The method of claim 20, wherein said earth formation 
includes at least one of rock, brine, hydrocarbon fxuid 
and clay. 

22. The method of claim 21, wherein said parameter which 
15 is determined forms the volume fraction of one of the 

hydrocarbon fluid and the brine. 
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